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Abstract 

We compute non-perturbative renormalization constants of fermio- 
nic bilinears for the chirally improved lattice fermions in the quenched 
approximation of QCD. We address finite size effects and the influence 
of Gribov copies. Our results are presented in the RT and MS schemes 
as well as in RGI form and we discuss relations between the renormal- 
ization constants implied by chiral symmetry. 

After publication we corrected the numerator of the first coefficient of 
a 3 s in (JMD from 3696847 to 3890527, which yields a 0.2% higher value 
of the conversion coefficient at \i = 2 GeV. 

PACS: 11.15.Ha, 12.38. Gc 
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1 Motivation 



Lattice Dirac operators obeying the Ginsparg- Wilson constraint [TJ (GWC) 
provide a lattice formulation of chiral symmetry closest to the continuum 
form; its violation is local and O(a). Exact GW-fermion fields are protected 
by their chiral symmetry [2] and renormalization of operators constructed 
from them is simpler than for e.g. Wilson fermions. Chiral symmetry implies 
several relations between renormalization constants. Checking these relations 
provides an important check of how well chiral symmetry is implemented in 
a particular lattice formulation of QCD [3]. 

To be more explicit, let us consider the local, flavor non-singlet quark 
field bilinear operator 

O r = uTd. (1) 

Here V denotes a Clifford algebra matrix. According to their Lorentz sym- 
metry we denote the five types of Y by S, V, A, T and P corresponding to 
scalar, vector, axial vector, tensor and pseudoscalar (S~ 1, V~ 7^, A~ 7^75, 
T~ Ih^i, Jv], P~ 75)- Chiral symmetry implies in particular Zs = Zp and 
Zy = Za- For conserved covariant currents, Ward identities give Zy = 1. 

So far only the overlap action [U [5] provides exact GW-fermions. There 
are, however, several lattice Dirac operators obeying the GWC approxi- 
mately (such as fixed point [6] and chirally improved (CI) fermions [7J) or 
in some limit [El Ej. Although technically more demanding than standard 
lattice fermion formulations they are still substantially less expensive (in 
terms of computer resources) than exact overlap fermions. The Bern-Graz- 
Regensburg (BGR) collaboration [TO] has been studying the low-lying hadron 
spectrum and other properties for fixed point [6] and CI fermions [TJ [H] . 

Results involving renormalized quantities like the pion decay constant, 
the chiral condensate or the quark masses, need renormalization constants 
in order to be connected to experiment, usually given in the framework of 
the MS (modified Minimal Subtraction) renormalization scheme. Here we 
present and discuss results for the renormalization factors of bilinear oper- 
ators of type ([T]) for the CI Dirac operator. A subset of our results, but 
without discussion of the full calculation, has been used earlier [T2] , 

The chirally improved Dirac operator has been introduced in [Til US] as 
an approximate solution of the Ginsparg- Wilson constraint. It is formulated 
as a truncated series of interaction terms with coefficients depending on the 
gauge coupling. Within the BGR collaboration it has been used to obtain 
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the ground state hadrons in the quenched approximation [TU] down to pion 
masses of 270 MeV, and it has been demonstrated to have good chiral prop- 
erties [H]. The renormalization constants for quark bilinears determined in 
this work permit a more detailed analysis of the chiral properties of the CI 
operator. 

In Sect. [2] we briefly discuss the method we use for calculating the renor- 
malization constants [T3] and its implementation. We then summarize the 
technical details such as the parameters of the simulations, the chiral limit, 
volume dependence and gauge fixing ambiguities in Sect. El before we present 
the lattice results in Sect. HI In order to compare with continuum notations 
like the MS scheme and to discuss the renormalization scale dependence, we 
recapitulate the conversion to other renormalization schemes in Sect. [5] and 
then discuss our results for the converted renormalization constants. We 
conclude and summarize in Sect. [61 



2 Method 

We want to compute renormalization constants non-perturbatively on the 
lattice. For this purpose we need a renormalization scheme which can be 
implemented in lattice Monte Carlo simulations as well as in continuum per- 
turbation theory. The latter property is necessary to enable the conversion of 
the lattice results to a more conventional scheme such as MS. A renormaliza- 
tion scheme satisfying these requirements is the so-called RI (regularization 
independent) scheme suggested by Martinelli et al. [13] . 

In this scheme one studies expectation values of the bilinear quark oper- 
ators between quark fields at a specific momentum value p 2 = fi 2 , 

(P I °r | p) \ p2=ti2 (2) 

and matches them to the corresponding tree-level matrix element (p | Or | 
p)q. This procedure is expected to work in a window 

A| CD « ^ « V« 2 (3) 

where discretization effects can be neglected, because the renormalization 
scale /i is small compared with the lattice cut-off 1/a, and (few-loop) contin- 
uum perturbation theory can be used to connect different schemes, because 
fi is much larger than the QCD scale parameter Aqcd- For comparing with 
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the MS scheme a typical value is // = 2 GeV. In most calculations one has 
a \i ~ 1 (or even somewhat larger) and the upper limit is not strictly obeyed. 
On the other hand, the limit also depends on the scaling properties of the 
actions involved. 

Since ([2]) is gauge- variant, one has to work in a fixed gauge and must 
compare the gauge dependent lattice matrix elements with the continuum 
results in the same gauge. Landau gauge fixing is a suitable choice, but 
one has to keep in mind that the Gribov copies uncertainty could spoil the 
comparison. In the lattice calculations one finds little, if any signal of such 
an effect (see our discussion below in Sect. [3] as well as [16j [TTJ HE]). 

Let us briefly summarize the method following [T5] in the modification of 
When multiplying (J2J) with (p | Op | p)$ 1 and taking the trace (note 



that the object in Eq. (j2]) is a matrix in color and Dirac space) one obtains 
for the renormalization condition 

Zr ^ tr[ (p | O r \ p) (p \ O r \ p)^ 1 ] I = 1 . (4) 

The matrix element 

(p\O r \p) = ^Ar(p) (5) 

is proportional to the amputated Green function 

Ar(p) = S-\p) G r (p) S^ip) , (6) 

and Z q is the quark field renormalization constant to be discussed below. 
The Green function Gr{p) is determined as the expectation value 

Gt(pU = \? E z- ip{x - y) E °M Mv)) ■ (7) 

' x,y \ z I 

The indices a, (3 run over color and Dirac indices and V denotes the lattice 
volume. The quark propagator is 

S a p(x,y) = (u a {x)up{y)) = (d a (x) dp{y)) (8) 

(assuming that u and d have equal masses and using the Landau gauge for 
the expectation value). 

So we have to compute Gr{p) and S(p). This is done in the following 
way. For the quark propagator S n evaluated on a single gauge configuration 
n we define 

Sn(x\p) = Y,e ipy S n (x,y). (9) 
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Taking into account 75-hermiticity of the propagator we may, for quark bilin- 
ear operators Or as defined in Eq. (QQ), rewrite Gr(p) in terms of the quantities 

w, 

Gr( P ) = i e~ ipix ~ y) (u{x)u{z)Td{z)d{y)) (10) 



1 N 



EE 75 s n {z\p) ] 75 rs n (z\ P ), (11) 



VN n=1 , 

averaging over N gauge configurations. Similarly we find the quark propa- 
gator in momentum space 

1 N 

s (p) « ^EE^W). (12) 

V ly n =l x 

S n (y\p) is computed by solving the lattice Dirac equation (Dei denotes the 
chirally improved Dirac operator) 

J2D CI (z,y)S n (y\p) = (13) 
y 

with a momentum source (cf. [19|). This has the disadvantage that one has 
to determine the quark propagators for several momentum sources, whereas 
in the original method [15] one uses point sources (i.e. taking into account 
just z = instead of summing over all z) and projects the quark sink to the 
desired momentum values. However, using momentum sources has the big 
advantage of a significantly better signal. 

The quark field renormalization constant is obtained by comparing the 
quark propagator to the free (lattice) propagator. Using the so-called RT 
scheme we take 



7 i _ 1 fr (c-i M IR{p)-ii^(pY 



(14) 



Here R and CL^ £1X6 the scalar and vector terms appearing in the free CI Dirac 
operator, which in momentum space reads Dci{p) = i^^a^ip) + R(p). Using 
Fourier transformation one can compute R and a M from the definition and 
the parameters in Dqi [3 [33] ■ They are normalized such that one finds 

a^ip) = ip^ + 0{ap) 2 and R(p) = 0{apf . (15) 
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Landau gauge fixing d^A^ = is implemented as discussed in [2U1 HZ] by 
iteratively minimizing a functional of the link variables with stochastic over- 
relaxation [21]. As is well known this type of gauge fixing still allows for 
Gribov copies. This gives rise to an uncertainty which we address in Sect. 

It is easy to check that at tree level one finds (p \ Or \ p)o = T. Putting 
things together we obtain the final formula for Z-p in the RF scheme 

r tr[A r (p)T-i] 

For T = 7^,7^75 averaging over the index v under the trace is implied. The 
relation between different schemes is discussed in Sect. El 



(16) 

j2 — * i.2 



3 Technicalities 

3.1 Parameters 

The gauge configurations were generated in the quenched approximation with 
the Liischer-Weisz action [22] at values of the gauge coupling (3 = 7.90, 8.35 
and 8.70 corresponding to lattice spacings of a = 0.148 fm = 0.750 GeV -1 , 
a = 0.102 fm = 0.517 GeV" 1 and a = 0.078 fm = 0.395 GeV" 1 , respectively 
[25] . The scale has been determined from the Sommer parameter r . 

The final numbers we quote were computed on lattices of size 16 3 x 32. 
For studying possible volume and gauge fixing dependence we used additional 
ensembles on 8 3 x 24 at /3 = 7.90. More details on the simulations can be 
found in [TO] . 

For each lattice size and gauge coupling we determined the quark propa- 
gators S(x\p) on different (gauge fixed) gauge configurations and for typically 
16 different momentum sources. We used 11 different values for the quark 
masses, ranging from am = 0.01 up to 0.20 covering a range of pion masses 
from 250 MeV up to 1-2 GeV (depending on 0). The boundary conditions 
for the fermions were chosen periodic in space and anti-periodic in the time 
direction. The momentum values were taken to lie roughly along the diago- 
nal of the Brillouin zone ranging from the origin up to p = (5, 5, 6, 10) 2n/L 
with values of a p up to 4.18. 

For each mass and each operator Or we evaluated Z^ 1 ' from f lT6|) . On 
the small lattice we used 10 gauge configurations and on the larger ones we 
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had 5. We calculated the statistical error with the help of the statistical 
bootstrap method. 



3.2 Extrapolation to the chiral limit 

As we want to obtain a mass-independent scheme we extrapolate all results 
for Z^ v (fi 2 , m ) linearly to the chiral limit m — >■ 0. In Fig. [T]this extrapolation 
is shown for three representative values of /i 2 . Whereas the linear fit is 
appropriate for most cases, we do find a substantially different behavior for 
the pseudoscalar density. Similar observations have been discussed in 



Since the pseudoscalar density couples to the Goldstone boson channel 
there are Oi\jm) contributions to Gp in the chiral limit. We can incorporate 
this fact by expanding the inverse renormalization constant [24] in terms of 
the mass and adding the pole term; in [23] only a pole term and a constant 
term were used, but one can also add terms 0(m) and higher, which vanish 
in the chiral limit (cf. [26]). We use 



1 ^ A(ji 



2\ 



Zp(j2 2 ,m) m 



+ B(fi 2 ) + C(^)m. (17) 



By subtracting the pole term first and then performing the chiral limit we 
can define a "subtracted" renormalization constant 

1 / 1 A(fi 2 



, r< , , - lim — — - = lim (B(fi 2 ) + C(fi 2 )m) = B(fx 2 ) , 

(18) 

such that the desired value is 

W) - ^ • (W) 

The fit to Zp on the large lattices, as shown in Fig. [TJ is very satisfactory. 
For the small lattice, which we use only to study the finite size dependence, 
the fit is not as good for the smaller values of /i 2 (cf. Subsection I3.4D . 

For large \i the operator product expansion guarantees the suppression 
of the pole contributions [15]. This behavior is clearly exhibited in the plot 
for Zp when comparing the results for smaller /i 2 to larger values of fi 2 . In 
terms of the expansion coefficients this means that A(fi 2 ) approaches fast 
enough for higher fi 2 . The appropriate order of limits would be first \i — > oo 
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Figure 1: Results for = 7.90 and volume 16 3 x 32: The values for Z' q 
(the fermion field renormalization constant), Z$, Zy, Zp, Za, Zp are plotted 
versus the quark mass at a 2 /! 2 = 0.163 (circles), 1.166 (triangles), 4.558 
(squares) together with extrapolations to the chiral limit. For Zp the non- 
linear fit function ( fl7l) is included (dashed line). The points at m = 
resulting from the corresponding subtracted extrapolation Zj, ub (/i 2 ) are also 
shown; the two values for the lower a 2 /i 2 overlap. 
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(a — > 0), then m — >■ 0. Here, due to the requirements of a lattice calculation, 
we first perform m — >■ 0. 

The axial vector also couples to the Goldstone boson. However, the cou- 
pling is proportional to the momentum transfer [27], which vanishes for our 
kinematics. 

For the further analysis we use the values in the chiral limit. 



3.3 Gauge fixing dependence 

Let us now analyze the effects due to gauge fixing ambiguities. As discussed 
above we fix our gauge configurations to Landau gauge using overrelaxation. 
If we apply a random gauge transformation to the original gauge configura- 
tion and then fix the gauge again, the procedure often leads to a different 
(gauge-equivalent) configuration. When comparing the results for the renor- 
malization constants from the two gauge copies we find small differences. 
This non-uniqueness is the Gribov uncertainty. We study the fluctuations 
induced by Gribov copies by comparing results for Z^ v (fi 2 ,m) derived with 
and without additional random gauge transformations before applying the 
gauge fixing procedure. In Fig. [2] we plot the ratio 

(where the superscript denotes the results from the two "different" configu- 
ration ensembles) as obtained for 8 3 x 24 lattices at /3 — 7.90 with Z % v taken 
in the chiral limit. For the values of fi 2 of interest (around 4 GeV 2 ) we find 
these effects to be substantially smaller than 1%. This result agrees with 
earlier studies [16l [18] . 



3.4 Finite size effects 

We use two different lattice sizes in order to study the volume dependence 
of the renormalization constants. In Fig. [3] we compare the values (extrapo- 
lated to the chiral limit as discussed above) for the RT scheme at (3 = 7.90. 
The deviations are generally small, increasing towards smaller momenta as 
expected for finite size effects. However, above /i 2 ^ 4 GeV 2 the difference 
is generally less than 1% except for Zp nh . A possible explanation of this 
exceptional behavior might be that the coupling to the light pions enhances 
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Figure 2: Values for the ratios Qs, Qv, Qt, Qa, defined in Eq. (I20l . vs. 
Ai 2 [GeV 2 ] for the 8 3 x 24 lattice at p = 7.90. 



the finite size effects for Zp nh . As mentioned in the discussion on the chiral 
limit, the extrapolation for Zp on the small lattices, in particular for small 
values of /x 2 , is not as good as on the larger lattices. This explains the larger 
error bars for lower /i 2 in this case. 

We restrict ourselves to the results from the largest lattice of size 16 3 x 32 
from now on. 



4 Results for RT 

In the appendix we collect our results for the chirally extrapolated values of 
Z^ 1 in tabular form. Figs. 0H6] summarize these data. We discuss in this 
section only the RI' values (triangles in Figs. 
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Figure 3: Values for Z' q , Z s , Z v , Z T , Z A , Z% ah vs. fi 2 at = 7.90. The 
triangles mark the values for the 16 3 x 32 lattice and the circles the values 
for the 8 3 x 24 lattice. Zf, ub is the value determined according to Eq. ( 1191) . 
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Figure 4: Values for Z' q , Z s , Z v , Z T , Z A , Z% ub vs. /i 2 [GeV 2 ] for the 16 3 x 32 
lattice at /3 — 7.90. RF values (triangles), MS values (squares) and RGI 
values (circles) are displayed. The dotted vertical line indicates the value 
/i — 2 GeV where the MS value has been obtained. Note that for the vector 
and axial vector the MS and the RGI values coincide. 
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Figure 5: Values for Z' q , Z s , Z v , Z T , Z A , Zf, ub vs. /i 2 [GeV 2 ] for the 16 3 x 32 
lattice at /3 — 8.35. RF values (triangles), MS values (squares) and RGI 
values (circles) are displayed. The dotted vertical line indicates the value 
/i — 2 GeV where the MS value has been obtained. Note that for the vector 
and axial vector the MS and the RGI values coincide. 
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Figure 6: Values for Z' q , Z s , Z v , Z T , Z A , Zf, ub vs. /i 2 [GeV 2 ] for the 16 3 x 32 
lattice at /3 — 8.70. RF values (triangles), MS values (squares) and RGI 
values (circles) are displayed. The dotted vertical line indicates the value 
/i — 2 GeV where the MS value has been obtained. Note that for the vector 
and axial vector the MS and the RGI values coincide. 
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Vector and axial vector have finite renormalization factors and should 
be scale independent. Indeed the changes with [i 2 (and also with the lattice 
size) are distinctly smaller than for the scalar or pseudoscalar renormalization 
factors. We also compared the variation of the renormalization factor for 
individual components of and and the components of T^ u . For fi 2 > 
4 GeV 2 the deviation from the mean value is smaller than 1.5 % for all gauge 
couplings studied, for fi 2 > 10 GeV 2 it is less than 0.3%. We therefore plot 
and list the average over the results for these components according to the 
remark following (|T6|) . 

Our currents are point-like and not the conserved currents. However, the 
values for the renormalization constants for vector and axial vector are both 
very close to 1, a behavior expected for conserved currents. 

The scalar and pseudoscalar densities are logarithmically divergent and 
have scale dependent renormalization. Their ratio, however, should be scale 
independent. Fig. [7] displays both ratios, Z^/Zy and Zp uh /Zs- Both ought 
to be 1 for chirally symmetric fermion actions. Indeed, the observed scale 
dependence of Zp nh and Z$ cancels in the ratio. Both ratios are surpris- 
ingly close to 1 with very little scale dependence, indicating the excellent 
approximation to chiral symmetry of the considered CI fermion action. 

Only for the two ensembles with smaller lattice spacing ((3 = 8.35, /3 = 
8.7) we have a reasonable chance to satisfy the limits fl3]). However, the 
relatively flat shape of the scale independent Za and Zy and the analysis of 
the \i dependence in the subsequent section indicate that we may still rely on 
the results around \i ~ 2 GeV. In fact, due to earlier results [TO] we expect 
that the Dirac operator used here has improved scaling properties even for 
larger lattice spacings. 



5 Conversion to RGI and MS 
5.1 Different renormalization schemes 

The connection between different renormalization schemes is established us- 
ing continuum perturbation theory. In the two schemes MS and RF the 
renormalization factors are related by 

7MS7,.2\ 

R ^ - WW) ■ < 21 » 
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Figure 7: Ratios Za/Zv and Zf? b /Zs for the data in the chiral limit (lattice 
size 16 3 x 32). 
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Both, Zp F (/i 2 ) and Zf^ s (/x 2 ) are in general divergent in the continuum limit, 
the ratio stays finite, though. 

Note that the RI' scheme which we are employing differs from the RI 
scheme only by the definition of the quark field renormalization constant. 
The ratio of the latter in the two schemes 



7RI 7RI 



(22) 



Z ' q ^ zr 

agrees with the conversion factor for the vector and axial vector renormaliza- 
tion constants from RI' to MS since there is no further renormalization due 
to the protection by Ward identities. In the Landau gauge and 3-loop order 
one has [2S1 EH] 

(23) 



Z^__Zf_ 67 /a s \ 2 / 52321 




Here ( n is the Riemann zeta function evaluated at n. This equation refers to 
the average over the different directional components of the currents (under 
the trace in (I16p ). For the scalar and pseudoscalar one finds 

ZW_Z^ 1 < 16 a. , f4291 152Cs r , ,, 

~ Zf ~ + 3 4tt + ^ 18 3 M^J ul ' 




+ ^ 3890527 _ 224993C 3 + 2960^ j'a.y + ^ 4) 

In order to illustrate the orders of magnitude involved we give a numerical 
value at the coupling a s (/x = 2 GeV) = 0.203576: 

« 1.16165 . (25) 

The correction for the different definitions of the quark field renormalization 
is rather small (~ 0.5%), whereas even the oc s contribution to the scheme 
conversion is large (« 2.9 %) in comparison. 

Concerning the tensors, the conversion factor from RI' to MS has been 
computed in [30]. In Landau gauge one finds 



r/MS 

Zjrp 
"7RF 




(26) 



100400 
27 
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We compute the coupling a s from the 3-loop expression in the MS scheme 
a s (q 2 ) 1 1 loglog(g 2 ) 



(27) 



J -an K} » (ft lo s 2 M<? 2 ) - ft lo s lQ g(^ 2 ) + & MS /?o - ft 



4vr /3 log(g 2 ) (31 logV 

1 

'/3 5 log 3 (g 

where we have introduced the shorthand notation q 2 = /i 2 /(Aqc D ) 2 . For 

Aqc D in the quenched approximation we use the value Aq^ d = 0.238 ± 
0.019 GeV calculated in [31]. 

In Figs. H]|6]we plot our results converted to the MS scheme. Except for 
the vector currents we expect dependence on both, /i 2 and lattice spacing. 

In Tables [TH the MS values are given at n — 2 GeV determined by linear 
interpolation between the closest \i values of the data. 



5.2 RGI values 

The renormalization factors depend on the renormalization scale \x (except for 
the vector and axial currents). This dependence is governed by the anomalous 
dimension 7 of the operator: 

7(a ,).-,i-.„Z« = | 7i (^) ,+1 . (28) 

Integrating this differential equation we can define a scale-independent quan- 
tity Z RG1 (RGI = Renormalization Group Invariant) by [32J 

(29) 

where the /3 function is given by 

« a , ) = ^ a , W = _4.g ft (^) ,+2 . (30) 

Since we know the renormalization group functions /3 and 7 only to a certain 
order in perturbation theory we cannot expect exact scale invariance when 
we apply the above definitions to our data. However, we can check how well 
the expectation of /i-independence is met [19J. 
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In Figs. HE we also show the RGI values as obtained by converting at 
each value of [i 2 (indicated by circles). For perfect 3- loop scaling these values 
should be constant; they are not, but they have indeed a plateau-like behavior 
for 2 GeV < fi < 4 GeV with less than 5% variation (10% for Zf G1 ). 

In Tables [T]|3] we quote the RGI values at afi = 1.08 (the actual values of 
a\i closest to 1). 

6 Summary and conclusion 

We have determined the renormalization constants of bilinear quark oper- 
ators for the chirally improved lattice Dirac operator by non-perturbative 
methods. These allow us to relate certain lattice observables to those in 
continuum renormalization schemes. We observe reasonable agreement with 
3-loop renormalization group behavior, which improves as the lattice spac- 
ing becomes smaller. Finite size effects and problems from Gribov copies 
were found to be negligible. The Goldstone boson pole in the renormal- 
ization factor Z P of the pseudoscalar density could be subtracted leaving 
behind a well-behaved Zp uh . The renormalization constants for the vector 
and the axial vector currents are close to 1. The ratios Za/Zv ~ 1.03 and 
Zp nh /Zs ~ 0.95 for the smallest lattice spacing indicate a very good approx- 
imation of chiral symmetry. 
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A Tables of results 



/i 2 [GeV 2 ] 


Z s 


Z v 




Z A 


Zp uh 


0.017 


0.7(1) 


0.65(8) 


1.0(1) 


1.1(1) 


0.98(2) 


0.291 


0.878(6) 


0.869(4) 


1.209(6) 


1.180(4) 


0.878(9) 


0.702 


0.84(1) 


0.887(2) 
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0.98398(2) 
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0.815(9) 


MS 
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0.9944(3) 


1.0087(4) 


1.0281(5) 



Table 1: Values for the renormalization constants in the RT scheme for the 
16 3 x 32 lattice at j3 — 7.90 as resulting from the extrapolation torn = 0. The 
RGI value is taken at /i 2 = 2.073 GeV 2 , i.e. a// = 1.07992. The MS (2 GeV) 
value is computed by interpolation. 
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1.012(1) 


0.987(4) 



Table 2: Values for the renormalization constants in the RT scheme for the 
16 3 x 32 lattice at /3 — 8.35 as resulting from the extrapolation to m — 0. The 
RGI value is taken at /i 2 = 4.365 GeV 2 , i.e. afi = 1.07992. The MS (2 GeV) 
value is computed by interpolation. 
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Table 3: Values for the renormalization constants in the RT scheme for the 
16 3 x 32 lattice at {3 — 8.70 as resulting from the extrapolation to m = 0. The 
RGI value is taken at li 2 = 7.464 GeV 2 , i.e. olli = 1.07992. The MS (2 GeV) 
value is computed by interpolation. 
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